Transport properties of nonhomogeneous segregated composites 
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In conductor-insulator composites in which the conducting particles are dispersed in an insulating 
continuous matrix the electrical connectedness is established by interparticle quantum tunneling. 
A recent formulation of the transport problem in this kind of composites treats each conducting 
particle as electrically connected to all others via tunneling conductances to form a global tun- 
neling network. Here, we extend this approach to nonhomogeneous composites with a segregated 
distribution of the conducting phase. We consider a model of segregation in which large random 
insulating spherical inclusions forbid small conducting particles to occupy homogeneously the vol- 
ume of the composite, and allow tunneling between all pairs of the conducting objects. By solving 
numerically the corresponding tunneling resistor network, we show that the composite conductivity 
a is enhanced by segregation and that it may remain relatively large also for very small values of 
the conducting filler concentration. We interpret this behavior by a segregation-induced reduction 
of the interparticle distances, which is confirmed by a critical path approximation applied to the 
segregated network. Furthermore, we identify an approximate but accurate scaling relation permit- 
ting to express the conductivity of a segregated systems in terms of the interparticle distances of a 
corresponding homogeneous system, and which provides an explicit formula for a which we apply 
to experimental data on segregated Ru02-cermet composites. 

PACS numbers: 64.60. ah, 73.40.Gk, 72.80.Tm, 72.20. Fr 



I. INTRODUCTION 

The transport properties of two-phase heterogeneous 
materials are strongly related to the structure of the com- 
posite i.e. the volume fraction <j) of the conductive fillers, 
their size and shape, and their dispersion into the insu- 
lating medium.— £ Controlling the conductivity a of dis- 
ordered composite materials by tuning these structural 
parameters is of fundamental importance for several ap- 
plications in the fields of micro and nanoelectronics such 
as, to name a few, electromagnetic interference shicld- 
ers, resectable fuses, strain and chemical sensors, flexible 
conductors, and anti-static compounds^ 

In several classes of composite materials in which 
the conductive particles are dispersed into a insulating 
medium, like e.g. polymer-based composites or metal- 
glass cermets, the prominent transport mechanism is 
quantum tunneling. In this case of tunneling, the con- 
ductance between two particles decays exponentially with 
the inter-particle distance over a characteristic tunnel- 
ing length £ which is of the order of a fraction to a few 
nanometers depending on the material characteristics. 

Typically, the tunneling mechanism is approximated 
by treating the fillers as core-shell objects, where the 
impenetrable hard core represents the physical particle 
and the thickness of the concentric penetrable shell is 
identified with ~ £. 4 ~ 7 In this way, the overall behav- 
ior of the composite conductivity a as a function of <fi 
is commonly interpreted in the framework of percolation 
theory*^ which considers two given conducting particles 
as either electrically connected or disconnected if their 
mutual distance is respectively lower or larger than the 
cut-off length imposed by the shell thickness. Under this 
assumption, the system undergoes a conductor-insulator 
transition at a critical concentration <p c of conducting 



phase, which is associated with the formation of a glob- 
ally connected cluster of electrically linked filler particles 
which spans the entire sample. Below <f> c there is no such 
sample-spanning cluster and a is zero, while for (f> > <p c 
the conductivity follows a power law behavior of the form 
a oc ((f) — 4>cY ■> where t is a critical exponent. 

Despite of its simplicity and of the unquestionable in- 
sights that it provides on the transport problem,— ~— in- 
troducing the concept of a cut-off through which elec- 
trical connections are established, is nevertheless a too 
crude approximation which alters the real nature of the 
interparticle electrical connectedness. This is indeed 
characterized by the fact that, due to the tunneling mech- 
anism, two given particles are always electrically con- 
nected regardless of their mutual distance^ even if the 
strength of the connection (i.e., the interparticle conduc- 
tance) decays exponentially with such distance. 

Recently, we have reformulated the transport problem 
in conductor-insulator composites by allowing each con- 
ducting particle to be connected to all others via tunnel- 
ing processes?^ - — and so without imposing the restric- 
tive hypothesis on which the usual core-shell model is 
based. By explicitly taking into account the composite 
morphology or microstructure and the conducting parti- 
cles shape and dimensions, this global tunneling network 
(GTN) model is able to describe the overall conductiv- 
ity dependence upon the particle concentration for many 
different classes of composites, ranging from granular- 
like systems to colloidal dispersions in a continuum in- 
sulating matrixJ^ When applied to the many published 
data on colloidal nanotube, nanofiber, nanosheets, and 
nanosphere composites, this formalism has permitted to 
extract important microscopic properties, such as the 
tunneling decay length £, directly from the experimen- 
tal conductivity vs concentration curvesi 15 i 16 
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FIG. 1: (Color online) Two-dimensional representation of (a) an homogeneous and (b) segregated dispersion of the conducting 
particles in the continuum. The insulating (conducting) inclusions are represented by open (filled) circles of diameter D2 (D 1). 
The insulating fillers can penetrate each other, while the conducting particles are impenetrable with respect to themselves and 
to the insulating ones. The grey region represents the space available for placing the centers of the conducting spheres, (b) 
Example of a distribution of insulating and conducting particles with D2/D1 = 4. 



In this paper we apply the GTN model to composites 
where the dispersion of the conducting fillers in the in- 
sulating continuum is not homogeneous. Specifically, we 
consider systems in which the insulating phase forbids 
the conducting fillers to occupy large (compared to the 
particle size) volumes inside the material, thereby lead- 
ing to a segregated spatial distribution of the conduct- 
ing phasei 6 ' 18 i 19 In real composites, like e.g. Ru02-bascd 
cermets^ this is achieved when the size of the insulating 
grains is consistently larger then the one of the conduct- 
ing fillers, and thermal treatments, inducing softening 
and sintering of the insulating phase without large-scale 
mixing, lead to a segregated distribution of the conduct- 
ing particles in a continuous insulating matrix. Besides 
the already cited Ru02-based cermets, another impor- 
tant class of segregated conductive composites is that of 
polymer-based ones ; 21 i 22 where an inorganic or carbona- 
ceous conductive filler is mixed with significantly larger 
polymer particles and the resulting compound is molded. 

The study of the transport properties in segregated 
composites is important for many technological appli- 
cations where low filler concentrations are demanded in 
order to have high conduction regimes combined with 
the unaltered mechanical properties of the host insulat- 
ing medium, or to reduce the quantity of the conduc- 
tive phase when its cost is high. So far, the problem 
of conduction in segregated systems has been limited 
to the evaluation of the critical concentration <j) c within 
the percolation framework, and both lattic o 18 i 21 ' 23 and 
continuum&i^ models have evidenced that 4> c is (usually, 
see RefflTjt) lowered by segregation. In the following we 
shall go beyond the percolation approach by consider- 
ing the GTN scenario for the segregation problem and 
by solving numerically the tunneling network equations. 
We find that the composite conductivity a for fixed filler 



concentration can be strongly enhanced by the segrega- 
tion, in accord with the observed trends. Furthermore, 
we show that the calculated filler dependencies of a can 
be reproduced to a great accuracy by the critical path 
approximation^ which we find to follow a simple scal- 
ing law permitting us to provide analytical formulas for 
a((f>). When applied to experimental data of real segre- 
gated composites, our formulas can be used to estimate 
the degree of segregation of the composite and the value 
of the tunneling length £. 

The structure of the paper is as follow: in Sec. [TT] we 
present our GTN model for segregated composites and in 
Sec. IIIII we calculate numerically the composite conduc- 
tivities. In Sec. IIVI we present our results on the critical 
tunneling distance which will be used to approximate the 
numerical results of Sec. IIIII and to provide explicit for- 
mulas for the conductivity. These are applied in Sec. [V] 
to some previously published data of segregated compos- 
ites to extract the tunneling distance. Section fVTl is left 
to the conclusions. 



II. MODEL 

We model the conductor-insulator composite as de- 
scribed in Fig. [IJa) and (b) for the case of an homoge- 
neous dispersions of conducting particles in the contin- 
uum and for a segregated distribution, respectively. In 
Fig- DIb) the spherical particles of diameter D2 represent 
the insulating inclusions (e.g., the glassy frit particles in 
Ru02 cermets) while the conducting particles are mod- 
eled as hard spheres of diameter D\. The two kinds of 
particles are mutually impenetrable and, furthermore, we 
assume that the D2 spheres can penetrate each other in 
order to simulate for instance sintering and softening of 



3 



the insulating grains. Typically, as in Ru02 cermets, Di 
is as large as a few micrometers while Di ranges from tens 
to hundreds of nanometers, so that the regime D2 D\ 
is the one of practical interest. Keeping this in mind, 
we shall consider in the following also moderate values 
of D2 1 D\ to better appreciate the overall trends towards 
the D2/D1 3> 1 regime. 

The system is generated by first randomly placing 
the penetrable insulating spheres into a three dimen- 
sional cubic volume of side length L with a given num- 
ber density P2 — N2/L 3 , where N2 is the number of 
D2 particles and L is chosen to be at least one or- 
der of magnitude larger than D 2 (we assume periodic 
boundary conditions). Since the positions of the insu- 
lating spheres are uncorrelated, their fractional volume 
is <f>2 = 1 — exp(— V2P2) where i>2 = ttD^/G is the vol- 
ume of a single sphered After having placed the insu- 
lating spheres, N% conducting hard spheres are added to 
the system through random sequential addition (RSA), 
where random placing is accepted only if there is no 
overlap with the other D\ and D2 spheres. The RSA 
procedure is repeated until the desired volume fraction 
value 4>i = piV\, where V\ = irDf/6 and p\ = Ni/L 3 , is 
reached. However, since the conducting spheres cannot 
penetrate the insulating ones, the available volume for 
placing the centers of the Di particles is reduced by the 
factor— 

v* = (l-fo) {1+Dl/D2)3 (1) 

with respect to the total volume L 3 of the system. This 
defines an effective volume fraction 4>\ /v* , larger than 4>\ , 
for the conducting spheres. Hence, for small values of v* , 
the limit achievable through RSA can become much lower 
than that of the homogeneous limit at v* = 1, which is 
^max ^ 0.382^ and scales approximately as cj)f & *v*. For 
densities larger than the RSA limit we have considered 
cubically arranged initial configurations where particles 
overlapping with the insulating spheres were removed. 
Both the initial RSA and the cubic configurations were 
then relaxed via Monte Carlo runs, where random dis- 
placements of the D\ sphere centers were attempted and 
accepted only if they did not overlap with any of its neigh- 
bors and with the D2 particles. Equilibrium was consid- 
ered attained when the mean nearest-neighbor distances 
between the D\ particles did not change within statisti- 
cal errors upon further Monte Carlo displacements. An 
example of the so-obtained distribution of D\ and D2 
spheres is shown in Fig.QJb) for the case D2/-D1 = 4. 

In describing the overall conductivity arising from the 
system described above we go beyond the usual core-shell 
approximation, and employ the GTN model to the subset 
of D\ particles of the composite. Hence, we treat any two 
conducting spheres centered at r, and rj as electrically 
connected through tunneling processes, irrespectively of 
their mutual distance = |r.; — r 3 -|. By assuming that 
the particle size D\ and the temperature are large enough 
to neglect charging energy effects, then the inter-particle 
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FIG. 2: (Color online) Histograms for the distribution P(er) 
of the conductivity obtained from 500 realizations of the tun- 
neling network for D 2 /D 1 = 8, = 50 and 4>i = 0.1. The 
volume fraction values of the insulating spheres are <j>2 = 
(with L/Di = 25) and 4> 2 = 0.4, 0.6 (with LjD x = 80). The 
solid lines are fits to log-normal distribution functions. 



conductance is given by: 

gij = g exp^-^-j (2) 

where go is a constant "contact" conductance which we 
shall set equal to the unity, £ is the tunneling decay 
length, and 5ij = rij—Di is the minimal distance between 
the surfaces of two conducting spheres. For a system 
composed by Ni particles, the GTN model is then equiv- 
alent to a weighted random network with N\ nodes, each 
with coordination number Ni — 1. However, contrary 
to the usual models of weighted networks,— the weight 
of each link is not random but it is given by Eq. 
which depends on the particular arrangement of the con- 
ducting fillers in the composite. This characteristics of 
the model permits in principle to study on equal foot- 
ing composites with different statistical properties of the 
microstructure, and has been successfully applied to ho- 
mogeneous colloidal and granular composites*^— In the 
following we shall show that the GTN approach is also 
able to describe the conductivity of segregated systems, 
thus providing a theoretical framework for the study of 
inhomogeneous composites. 



III. CONDUCTIVITY 

In this section we present the results of our numerical 
calculations of the filler concentration dependence of the 
conductivity a for different segregated systems specified 
by <f>2 and D2 /D\ . Since in the GTN model each conduct- 
ing particle is connected to all others through Eq. @, the 
calculation of a would require the solution of a network 
with Ni(Ni — l)/2 resistors, which is a computationally 
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FIG. 3: (Color online) Calculated GTN conductivity as a function of the volume fraction 0i of the conducting spheres with 
diameter D\ for = 15 and 50, and for different values of the volume fraction 02 of the insulating spheres with diameter 

D2. (a) D2/D1 = 4, (b) D2/D1 — 8. The solid lines in (b) are fits of our formula. 



demanding, or even insurmountable, task for the Ni val- 
ues considered in our study. However, we can exploit 
the exponential decay of Eq. ((2J) by neglecting contri- 
butions from tunneling between particles sufficiently far 
away apart. Indeed, depending on the value of £ and of 
the filler concentration 0i, it is possible to identify an 
upper artificial cut-off S* such that the conductances be- 
tween particles at mutual distances Sij > S* can be safely 
removed, reducing drastically the number of connected 
particles in the network. In our calculations we have 
then chosen S* such that, depending on 0i, exp(— 28* /£) 
is from five to twenty orders of magnitude smaller than 
the overall network conductivity. 2 - 

Once the network has been reduced, we have eval- 
uated a by combining the numerical decimation algo- 
rithm of Ref. with a preconditioned conjugated gra- 
dient method. Specifically, we decimated iteratively the 
network starting from the nodes with the lowest coordi- 
nation number in order to eliminate dead ends and to 
compact the network. We continued the decimation pro- 
cedure until a single conductance was left, whose value 
coincided with the conductance of the original network. 
If the computational time for the node decimation was 
too large, as it was typically the case for segregated sys- 



tems with 4>\/v* large, we switched to the conjugate gra- 
dient method (sec e.g. Ref. I30I ) applied to the partially 
decimated network. We have applied this procedure to 
N r = 200 — 600 realizations of systems with N% con- 
ducting spheres ranging from N± = 250 (for L/D\ = 50 
and 0i = 10~ 3 ) to iVi = 322690 (for L/D x = 80 and 
0i = 0.33). 

In Fig. [5] we show the distributions P(a) of the a val- 
ues obtained from 500 realizations of the system with 

01 = 0.1, D r /i = 50, D 2 /D 1 = 8, and for 2 = (with 
L/Dx = 25), 0.4, and 0.6 (with L/D 1 = 80). All three 
sets of data follow approximately a log-normal distribu- 
tion (solid lines) which stems from the exponential decay 
of Eq. © (note that the distribution of the 02 = case 
is broader than the two others because of the smaller size 
of L). The distributions arc peaked at the average of the 
logarithm of a, with no significant drifts when the sys- 
tem size is increased for fixed 0i. It is clear from the 
figure that, although 0i is kept fixed, the mean value of 
the conductivity steadily increases as the volume fraction 

02 of the insulating spheres is enhanced. From Eq. (fTJ) 
this trend can be interpreted by noticing that as 02 in- 
creases, the available volume fraction v* decreases, so 
that the conducting particles occupy a reduced volume 
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compared to the homogeneous 4>2 = case (i.e., the effec- 
tive concentration 4>\/v* is larger). In turns this means 
that, as it can be inferred from Figs. \V[a) and (b), the 
mean interparticle distances Sij in Eq. ([5]) are reduced for 
4>2 =/= 0, thereby leading to an enhancement of the overall 
conductivity. 

This behavior is clearly shown in Fig. [3l where we plot 
tr (symbols) as a function of (j>\ for D2/-D1 = 4, Fig.^a), 
and D 2 /Di = 8, Fig. [3fb), and for several values of 
the insulating sphere densities <p2- The reduction in a 
for decreasing 4>i, which is a direct consequence of the 
fact that as 4>i is reduced the interparticle distances S^ 
get larger, can be strongly mitigated by the segregation 
which, through the reduction of v* , tends instead to de- 
crease Sij. We have therefore that, as shown in Fig. [31 
for fixed values of the tunneling length £, as the segre- 
gation is enhanced the threshold value of <\>\ required to 
achieve a given a decreases considerably. By combining 
this result with the observation that, in practice, the low- 
est measurable conductivity a m - la in real composites is set 
either by the experimental setup or by the intrinsic con- 
ductivity of the insulating phase, and that this defines a 
sort of "critical" threshold <fii c at which a ~ own,— we 
obtain that more segregated systems entail lower values 
of <f>\ c . This last observation is in agreement with the 
behavior seen in real segregated composites if we rein- 
terpret the percolation threshold values reported in the 
literature as our crossover concentration <p\ c . 

As we shall see in the next section, the interpretation 
that the segregation basically leads to a shortening of 
the tunneling distances can be established on a firmer 
ground by employing the critical path approximation to 
the tunneling network. This analysis will also provides 
us with useful explicit formulas for the overall composite 
conductivity. 



IV. CRITICAL PATH APPROXIMATION 

We show in this section that, as for homogeneous dis- 
persions of spherical, rod-like, and plate-like impenetra- 
ble particles ; 15 ' 16 the critical path approximation (CPA) 
can also reproduce to a high accuracy the conductivity 
behavior for the inhomogeneous dispersions considered 
here. According to the CPA,— the composite conductiv- 
ity a can be expressed approximately as 

a = cr cxp^-|(5 c ^ , (3) 

where S c = S c ((f>i, <f>2, D2/D1) is the largest among the 
5^ distances such that the sub-network defined by all 
distances Sij < S c forms a conducting cluster span- 
ning (or percolating) the entire sample. The critical 
distance S c defines the single bond critical conductance 
g c = cxp(— 2<5 c /£) which, once assigned to all the con- 
ductances Qij of the network, leads to Eq. ([3]) , where the 
prefactor ctq is the only remaining fitting parameter. 
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FIG. 4: (Color online) Histograms for the distribution P(5) 
of the minimum interparticle distance <5 required to have per- 
colation for 500 realizations of the same systems considered 
in Fig. [2] The solid lines are fits to gaussian distribution 
functions. 



Equation @ reduces the conductivity problem to a 
simpler, geometrical one, which amounts to find the ge- 
ometrical critical distance S c so that percolation is es- 
tablished. In practice, S c can be obtained by coating 
each conducting sphere by a concentric penetrable shell 
of thickness 5/2 and by considering two spheres as con- 
nected if their shells overlap. The critical distance S c is 
then the minimum value of S such that (for given val- 
ues of (pi, (f>2, and D2/D1) a cluster of connected spheres 
spans the sample. Hence, in contrast to the usual core- 
shell model approach to transport,— "—1^— here the shell 
thickness is not fixed a priori but it depends on the par- 
ticle concentration. 

Our numerical procedure to find S c goes as follows. For 
fixed volume fraction (pi of the conducting spheres, as 
well as for given fa and D2/D1 values, we first generate 
the system as explained in Sec. [Ill For each realization i 
of the system (i = 1, . . . , N r ) we chose an initial value Si 
comprised within the interval AS = <5 max — <5 m m, where 
Smhi = and <5 max is large enough so that a percolating 
cluster is surely established. Clustering is performed on 
the adjacency list which represents the vicinity network. 
Namely, we scan iteratively the list by checking if a given 
node already belongs to previously classified clusters. If 
this condition is not fulfilled, we form a new cluster by 
identifying which nodes are directly and indirectly con- 
nected to the selected node. This is done by labeling the 
first and, recursively, the next levels in the vicinity hier- 
archy. Finally, a percolating cluster is the one in which at 
least two of its nodes lie at the opposite faces of the sam- 
ple cube. The critical distance S l c for the i-th realization 
is then found by bisecting the interval AS until conver- 
gence is reached within a relative error of 10 -3 . Finally, 
an histogram representing the distribution P(5 C ) of the 
critical distance is obtained by repeating the procedure 
for all the N r realizations of the system. 
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Examples of the thus obtained P(8 C ) arc shown in 
Fig. H] for the same cases of Fig. [2] (i.e., (f>i = 0.1, 
D 2 /D 1 = 8, and 2 = 0, 0.4, and 0.6). They approx- 
imately follow normal distributions (solid lines) centered 
at critical distances which steadily decrease as one moves 
from the homogeneous case {<f>2 = 0) to the increasingly 
segregated regimes (02 = 0.4 and 0.6), thus confirming 
our previous conjecture that segregation implies a short- 
ening of the tunneling lengths. 

This trend is clearly seen in Fig. [5j where we plot the 
overall behavior of 5 C as a function of <pi and for differ- 
ent values of 4>2 and D2/D1. At very low filler volume 
fractions, all S c curves for ^ ^0 tend asymptotically to 
the critical distance of the homogeneous case (open cir- 

— 1/3 

cles) which in this limit behaves as S c oc (f> 1 , indicating 
that segregation is irrelevant for <f>i — > 0. As <\>\ increases, 
segregation acts by lowering the critical distance which, 
for large <f>2 and D2/ D\ values, can be even one order of 
magnitude smaller than that of the homogeneous case. 

The reduction of S c can be explained by using the ar- 
gument that the available volume for placing the con- 
ducting spheres is reduced by segregation. Indeed, as 
Fig. OJb) suggests, such a reduction has the net effect of 
increasing the local density of the D\ particles, and so of 
reducing the critical distance for percolation. This argu- 
ment neglects conductor-insulator interface effects and 
assumes that narrow bottlenecks are irrelevant for the 
establishment of the critical paths, which are conditions 
both fulfilled in the D 2 /D 1 > 1 case^- 

By following this line of reasoning, we argue that the 
only relevant variable for describing segregation is the 
available volume fraction v* , and that consequently the 
critical distance can be expressed as S c (4>i, v*). Further- 
more, by requiring that the relevant 4>i dependence of 
the Di particles is through the effective volume fraction 
4>\/v*, we express the critical distance as S c (4>i,v*) = 
af((f>i/v*), where f(x) is a general function and a is a 
proportionality constant. The value of a can be deter- 
mined by requiring that the critical distance for <f>i — > 
coincides with that of the homogeneous case, which leads 
us to the following scaling relation 

8 c {4> l ,v*)=v*- 1 /*f{j> 1 /v*). (4) 

As shown in the insets of Fig. [51 where v*^ 3 S c is plotted 
as a function of <f>x/v*, the data for different <j) 2 values 
basically collapse into a single curve already for D2/D1 > 
8. However, for D 2 /Di =4 the scaling relation @ is not 
very effective but, as shown in the Appendix, the scaling 
argument can be generalized in order to provide a better 
data collapse also for D 2 /Di < 8. 

From the scaling relation Q it follows that for v* = 1 
the function f(4>i) coincides with the critical distance in 
the homogeneous limit 5 c (<fri,l), thereby leading to 



5 c (c/> 1 ,v*) = v*-^ 3 S c (<f> 1 /v*,l), 



(5) 



which merely states that the critical distance in the segre- 
gated regime can be directly obtained from that of the ho- 
mogeneous case. In order to illustrate the consequences 
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FIG. 5: (Color online) Critical distance S c dependence on the 
volume fraction cj>\ of the conducting D\ spheres for several 
values <j>2 of the volume fraction of the insulating D2 spheres 
for (a) D 2 /D 1 = 4, (b) D 2 /D 1 = 8, and (c) D 2 /D 1 = 12. 
Insets: the same data plotted according to the scaling relation 
of Eq. ©. 



of Eq. §5§ for the conductivity a of segregated compos- 
ites, let us first verify that the CPA of Eq. ([3]) actually 
provides a valuable approximation of a . In Fig. [5] we re- 
plot the conductivity data of Fig. [3] as a function of the 
critical distance results of Fig. [5j Irrespectively of 4>2 and 
of D2/D1, we find that a nicely follows the linear relation 



ln(tr) = ln(<T ) 



(6) 
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FIG. 6: (Color online) Natural logarithm of the conductivity values a of Fig. [3] plotted as a function of the critical distance S c 
of Fig. [5] (a) D2/D1 = 4; (b) D2/D1 — 8. The insets shows the ln(cr) versus 5 C /D\ dependence in a larger scale. 



with slope 2/£, and that therefore the CPA is in excellent 
agreement with the full numerical solution of the tunnel- 
ing resistor network. Thus, from Eqs. and ([5]), this 
means that for D2/D1 ^> 1 (see the Appendix for the 
D2/D1 ~ 1 case) we can express a in terms of the crit- 
ical distance 6 c (4>i, 1) of the homogeneous limit, leaving 
only the prefactor <jq to be determined. By using some 
approximate formula for 5 c (4>i, 1), like for example those 
reported in Refs. ll5lll6ll35l . the full </>i dependence of a for 
given v* can be then expressed in analytical terms with 
high accuracy. This is illustrated in Fig. Elb) , where the 
solid lines have been obtained from 

r 2 ^(^/^,1) 1 w 

a = ( 7 cxp^-- ^ 1/3 j, (7) 

with v* from Eq. ([1]), Dq,/D\ = 8, and 6 c (x, 1) as given 
in Ref. [LIU^ As shown in the next section, the possi- 
bility of expressing the composite conductivity in terms 
of an analytic formula is a valuable tool to describe the 
filler dependences of real segregated composites and to 
interpreted the experimental data. 
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FIG. 7: (Color online) (a) Experimental conductivity a data 
of two series of Ru02-cermets as a function of the volume 
fraction <fri of the Ru02 particles (from Ref. l32h . The same 
<r values are plotted as a function of the critical distance S c 
for (a) D\ ~ 100 nm and (c) D\ ~ 40 nm by assuming two 
different values of the available volume fraction v* . The solid 
lines are fits to Eq. ((6}. 



V. APPLICATION TO EXPERIMENTS 

We show here that the result of the previous section 
can be used to analyze experimental a vs <pi data in order 
to extract estimates of the tunneling length £ value and 
of the degree of segregation in the composite. This is so 
because, if the experimental values of ln(<r) are plotted 
as a function of 6 C instead of </>i , and if our GTN picture 
applies, then they are expected to follow Eq. ([6]) whose 
slope — 2/£ directly gives the value of the tunneling factor 
£ independently of the specific ctq value. Furthermore, 
the available volume fraction v* appearing in Eq. ([S]) docs 



not depend, in principle, on the particular model chosen 
to represent segregation,— and so it may used as a fitting 
parameter which best reproduce the experimental data. 

In order to illustrate how the theory applies to real 
composites, we consider here the conductivity data of 
RuCvcermet samples which were already reported in 
Ref. [321 . In particular we consider two series of samples 
constituted by Ru02 conducting particles of mean sizes 
D\ « 40 nm and D\ « 100 nm, dispersed in a borosili- 
cate glass. The glassy grains prior to thermal processing 



8 



10 
9 



„ 6 



£■ 5 



-■- D 1 ~ 100 nm 
-•-D^^nm 



0.2 



0.4 



0.6 
* 

u 



1.0 



FIG. 8: (Color online) Values of the tunneling decay distance 
£ as a function of the available volume fraction v* extracted as 
explained in the text. The two horizontal dashed lines at 4.6 
nm and 2.1 nm represent the maximum and minimum values 
of £ which are compatible with both series of composites. 



(firing) had average size of about 3/mi, so that for both 
series of composites D2/D1 ^> 1. The two series of sam- 
ples were fired by following identical thermal cycles so 
that, in principle, they differ only in the mean size D\ of 
the conducting Ru02 particles. It should be noted how- 
ever that although the finer RuO"2 powders were given by 
nearly spherical and monodispersed particles, the coarser 
powders had more dispersed grain sizes with less regular 
shape. 

In Fig. [7{a) we plot the measured conductivity as i\ 
function of Ru02 volume fraction </>i for both series of 
composites. In Ref . l32l we interpreted these same data in 
the framework of percolation theory and fitted them with 
the power-law relation (</>i — </>i c )*. The resulting low per- 
colation threshold values, <j)\ c ~ 0.07 — 0.05, were found 
to be consistent with the segregated distribution of the 
RuC>2 conducting phase observed in the microstructure, 
while the large transport exponent values, t ~ 3 — 4, were 
concluded to arise from the nonuniversality of the critical 
behavior as predicted by the tunncling-percolation model 
of Re f. KpH ( see also later developments of this theory in 
Refs. I3hl34r ). Here, we offer an alternative interpretation 
of these data based on the GTN theory which, as ex- 
plained in this paper and in Rcfs.[L5T-[l7l is more justified 
on physical grounds than our previous percolation-based 
one of Ref [32. 

Let us start by re-plotting the conductivity data of 
Fig- Eta) in terms of the critical distance 5 c {<f)\,v*) by us- 
ing the scaling relation of Eq. ([5]) . For the functional form 
of 8 c (6i/v* 1 1) we use the fitting formula published in 
Ref. I15U 6 - and we treat v* as an adjustable parameter^ 
The resulting ln(cr) vs 5 C plots for two different values of 



v* are shown in Fig.^b) and Fig.JTJc) for the D\ w 100 
nm and D\ sa 40 nm samples, respectively. In both cases 
the data follow a better linear dependence for v* = 0.3 
(strong segregation) than for v* = 1 (homogeneity), and 
from the corresponding fits to Eq. ([6]) (solid lines) we 
obtain that lower values of v* imply also lower values 
of £ (i.e., the slopes are larger). This is most clearly 
seen in Fig. [5] where we report the v* dependence of the 
so-obtained £ values for both series of samples. Start- 
ing from the homogeneous limit at v* = 1 the tunneling 
length £ of the D\ « 100 nm series decreases by a fac- 
tor of five when the available volume fraction is lowered 
down to the minimum value v* = 0.2 for which a linear 
fit of ln(er) vs S c was possible, while the D\ « 40 nm 
case displays a weaker decrease in the same range of v* 
due to the smaller Ru02 grain size. By realizing that 
the tunneling decay length should be independent of the 
size D\ of the conducting particles, while the specific v* 
value could be different for the two series of composites, 
then £ must be comprised between the two horizontal 
dashed lines at 4.6 nm and 2.1 nm in Fig. [5J However, 
since the microstructure of both series of composites dis- 
plays a marked segregated dispersion^, then v* should 
be sensibly smaller than the unity, suggesting that the 
lower limit £ ~ 2 nm is a more reliable estimate for £. 
This value is fully comparable to those extracted from 
other conductor-insulator composites ) 15 ' 16 and, specifi- 
cally, agrees well with the results of microscopical inves- 
tigations of thick-film cermet resistors^ 



VI. CONCLUSIONS 

In this paper we have generalized the GTN model, 
where each conducting particle is connected to all oth- 
ers through tunneling processes, to describe composites 
whose microstructure is given by a segregated dispersion 
of the conducting phase. This particular class of non- 
homogeneous composites is characterized by large con- 
ductivity even for volume fraction values <f>\ of the con- 
ducting phase as low as a few percents. According to the 
percolation theory, this behavior is explained by the re- 
duced values of the percolation threshold induced by the 
segregated dispersion of the conducting particles. Here, 
we have shown that the <j>i dependence of the conduc- 
tivity in nonhomogeneous segregated composites can be 
understood without imposing any fixed cut-off in the mi- 
croscopic electrical connectivity (as it is done in percola- 
tion theory) and that the GTN formulation provides thus 
a natural and physically justified approach to the study 
of transport in disordered composites. 

Besides the full numerical solutions of the tunneling 
resistor network, we have also shown that the critical 
path approximation is valid for a wide range of <pi, and 
that it permits to formulate a scaling relation connecting 
the critical tunneling distance 8 C for a segregated sys- 
tems with that of a homogeneous composite. Finally, 
we have illustrated the practical importance of this seal- 
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FIG. 9: (Color online) Dependence of the correction parame- 
ter (a) a and (b) /3 on the volume fraction (j)2 of the insulating 
spheres and on D2/D1. Inset: the 8 C data for D2/D1 = 4 of 
Fig.[5jb) rescaled according to Eq. ()A2[) for the corresponding 
values of a and /3. 



ing by applying it to experimental conductivity data of 
Ru02-cermet segregated composites, which has permit- 
ted us to extract a realistic tunneling decay length £ and 
to estimate the degree of segregation in these materials. 
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Appendix A: improved scaling formula 

Here we briefly address the problem of generalizing the 
scaling relation of Eq. ([5]) in order to extend its validity 
beyond the D 2 /Di ^> 1 regime. To this end, let us first 
remind that v* denned in Eq. ([1]) gives the volume frac- 
tion available for placing the centers of the conducting 



particles with diameter D\. However, in order to de- 
fine the critical distance S c , the conducting spheres are 
treated as core-shell particles, where the hard-core of di- 
ameter Di is coated by a concentric penetrable shell of 
thickness S c /2. Since this penetrable shell may actually 
overlap the insulating spheres then, to what concerns the 
connectivity, these latter may be treated as having effec- 
tively a smaller diameter aD 2 < -D2, where < a < 1 is 
a correction factor which captures such effective reduc- 
tion. In this way the corrected available volume fraction 
reads: 



(1-0: 



^(a+D 1 /D 2 ) 3 



(Al) 



For a < 1 the net effect of the overlapping between shells 
and insulating spheres is then an effective increase of the 
available volume. This increase is expected to be unim- 
portant when D2/DI ^> 1, so that a ~ 1 in this regime, 
while when the two diameters D\ and D 2 are comparable 
the a correction has to be considered. By using Eq. (|A1|) , 
we generalize thus the scaling relation (0| as follows: 



£c0i,O = w* /3 /(0i/<), 



(A2) 



where in addition to the corrected available volume f* 
we have introduced the new exponent /3 to improve the 
scaling in the large (f>\ region. 

In Fig. [S] we plot (a) a and (b) (3 obtained from the 
minimization of \\5 c (<pi/v^, : 1) — w* _,9 (5 c || as functions of 
the volume fraction (f> 2 of the insulating spheres and for 
D 2 /D l = 4, 8, and 12. For D 2 /D 1 = 8 and 12 the co- 
efficients a and /3 are close to respectively a = 1 and 
(3 = 1/3 for all values of (j) 2 considered, so that, as ex- 
pected, Eq. ([5]) with v* as given in Eq. ((T|) provides a 
rather good scaling of the 8 C data. On the contrary, for 
D 2 /Di = 4 the coefficient a displays a stronger (f> 2 de- 
pendence and is sensibly smaller than a — 1 for large <j) 2 
values, while /3 is larger than 1/3. In particular, a < 1 
indicates that for this case the increase of the effective 
available volume u* is an important effect for the correct 
scaling of 8 C which, as shown in the inset of Fig. (Hla), is 
now almost perfect. 
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